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Abstract 



We study the effect of virtual Coulomb break-up, commonly known as 
the dipole polarizability, of the deuteron projectile on the astrophysical fu- 
sion reaction ^He((i,p)^He. We use the adiabatic approximation to estimate 
the potential shift due to the E\ transition to the continuum states in the 
deuteron, and compute the barrier penetrability in the WKB approximation. 
We find that the enhancement of the penetrability due to the deuteron break- 
up is too small to resolve the longstanding puzzle observed in laboratory 
measurements that the electron screening effect is surprisingly larger than 
theoretical prediction based on an atomic physics model. The effect of the 
^He break-up in the ^He(d,p)'^He reaction, as well as the ^Li break-up in the 
'^Li(p, a)^He reaction is also discussed. 
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The problem of electron screening effect on nuclear fusion reactions measured at a lab- 
oratory at very low incident energies has not yet been fully understood. The rise of the 
astrophysical ^'-factor for reactions such as ^He((i,p)''He and D{^B.e, p)''He as the incident 
energy goes down below about 50 keV has been attributed to the screening effect of the 
bound electrons in the target atom (or molecule), which shields the Coulomb potential 
between the colliding nuclei [1,2]. It has been found, however, that the amount of the 
enhancement of the S'-factor can be accounted for only when an unrealistically large elec- 
tron screening energy t/g is used in a calculation [3-7]. The value of 1/^ required to fit the 
data ranges between 0.88 and 14.5 times the adiabatic value (see Table 1 in Ref. [8] for a 
summary), where the screening energy is given by a difference of electron binding energies 
between the unified and isolated systems [1,9]. Since the adiabatic approximation should 
provide the upper limit of the screening energy [10-13], the mechanism of the low energy 
enhancement in the ^'-factor remains an open problem. A noteworthy recent paper is by 
Barker, who refitted the experimental data by including the screening correction as a free 
parameter and obtained smaller screening energies which are consistent with the adiabatic 
value in many systems [14]. However, for some cases, the optimum screening energy still 
exceeds the adiabatic value, and the fitted S factors are somewhat in disagreement with the 
experimental result of the Trojan-horse method [15-17], which is believed to provide the 
bare cross sections without the influence of bound electrons. Thus, the problem has not 
been resolved completely yet. 

Besides the electron screening effects, several small effects on astrophysical fusion re- 
actions have also been examined. These include vacuum polarization [18], relativity [18], 
brcmsstrahlung outside the barrier [18], atomic polarization [18], radiation correction during 
the tunneling [19], zero point fluctuation of nuclei in the atom and the molecule [8], and the 
effect of flnite beam width [8] . All of these effects have been found much smaller than the 
screening effect. 

In this paper, we consider more corrections to astrophysical fusion reaction. An impor- 
tant fact is that the classical turning point of interest is much larger than the nuclear size 
(for instance, it is 288 fni for the d+'^He reaction at ii^c.m.=10 keV), and effects which are 
relevant to the reaction have to be associated with the Coulomb interaction or otherwise 
very long ranged. The effects associated with the nuclear interaction will be washed out by a 
careful choice of effective nuclear potential between the projectile and the target, unless the 
energy dependence is very strong [19]. In this sense, the nuclear absorption under the barrier 
[20-22], which has been discussed in connection with the sharp rise of nuclear S'-factor for 
the ^^C -I- ^^C fusion reaction at low energies, for instance, is not helpful for the astrophysical 
reactions. One may also think about the non-local effects of the internuclear potential on 
the tunneling phenomena [23-25]. However, this effect does not seem significant cither, as 
one can cast from the result of microscopic cluster model calculation for the d+^Hc reaction 
at low energies [26] , where the exchange effect has been included both in the nuclear and in 
the Coulomb interactions. 

As another example of small effect on fusion, we here consider the effect of Coulomb 
break-up of colliding nuclei. At energies which we are interested in, the break-up channel 
is most likely kincmatically forbidden. However, the tunneling probability is still infiuenced 
through the virtual process [27], and it is important to estimate the size of its effect for 
a complete understanding of the reaction mechanism. This effect is also known as dipole 



2 



polarizability. To our knowledge, this effect has not yet been computed in the hterature, al- 
though a few calculations have existed based on the continuum-discretized-coupled-channels 
(CDCC) method for transfer reactions at much higher energies, which were performed in 
aiming at extracting the cross section of the astrophysical radiative capture reactions at zero 
incident energy [28,29]. We also notice that the enhancement of tunnehng probability due 
to the break-up coupling has been extensively discussed in the context of subbarrier fusion 
reaction of a halo nucleus [30-34]. 

We use a three-body model in order to estimate the effect of the virtual Coulomb excita- 
tion of projectile on the tunneling probabihty. Denoting the coordinate between the target 
and the center of mass of the projectile by R and the coordinate between the projectile 
fragments by r, the Coulomb interaction in this system reads 

V (R r) — ZiZxe^ ^ Z^Zt^? 

*^ ' |i?-Fm2r"/(mi -Fm2)| |i2 - mir/(mi + 7712) | ' 

where T^^ is the El operator given by 

ff = eEirFi^(f). (3) 

Here, cei is the El effective charge given by {m2Zi — miZ2)e / {mi +m2), where mi and 1712 
are the mass of the projectile fragments while Zi and Z2 are their charges {Zp = Zi + Z2 is 
the total charge of the projectile). For a head-on collision, the incident channel is a s-wave 
bound state (t)Q{r) of the projectile coupled to the relative angular momentum L = for the 
R coordinate. This channel couples to a p- wave state 0i(r) of the projectile via the coupling 
interaction (2). The relative angular momentum for the R coordinate has to be L = 1 in 
the excited state channel so that the total angular momentum is conserved. For simplicity, 
we have neglected the spin of the projectile fragments. The matrix element of the coupling 
potential between these channels is given by 



where B{E1) |(0o||7'^^||0i)P is the strength of the electric dipole transition of the 
projectile. 

For an exponential wave function for the bound state 00 together with the plane wave 
function for the scattering state 0i, a simple and compact expression for B{E1) | has been 
derived by Bertulani, Baur, and Hussein [35], which is given by 

dBjEl) T ^ Sh'el, ^{E^ - ef/^ 
dE, 7r//?2 ' 

where //12 = mim2/(mi -|- 777.2) is the reduced mass of the projectile system and e is the 
binding energy. This function has a peak at E^ — 8e/5, and the total dipole strength is 
given by [35] 
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Bm T= (6) 

In this work, for simplicity, we assume that the El strength is exhausted by a single state 
at E^. = 8e/5 with the strength given by Eq. (6). With this prescription, the problem is 
reduced to the two dimensional coupled-channels calculation with the coupling matrix given 

by 



where /i is the reduced mass for the R motion. 

In order to estimate the coupling effect, we use the adiabatic approximation and derive 
the adiabatic potential shift by diagonalizing the coupling matrix (7) at each R. Taking the 
smaller eigenvalue, the potential shift is given by 



AVUR) = ^ ^\ ' ' ■ (8) 

Note that this potential shift coincides with the adiabatic polarization potential which Alder 
et al. derived using the second order perturbation theory [36] (see also Ref. [37] for a 
derivation using the Feshbach formalism), in the hmit oi E^ F{R) and when one ignores 
the angular momentum transfer. We then compute the tunneling probability using the WKB 
formula for low energy. 



P{E) = exp 



2 ' dRj'^iVoiR) + AVUR) - E) 



(9) 



where Vq{R) — ZpZtg^ / R'ls the bare Coulomb interaction, and Rq and Ri are the inner and 
the outer turning points, respectively. Notice that the adiabatic approximation provides the 
upper limit of the potential penetrability [10-12]. Therefore, our results should be regarded 
as the upper limit of the virtual break-up effects in the astrophysical reactions, although the 
adiabatic approximation should work well at astrophysical energies. 

The effect of the target break-up can also be taken into account in a similar manner. In 
this case, one considers five channel states: i) the incident channel, ii) the projectile break-up 
channel, iii) the target break-up channel, iv) the mutual break-up channel with the relative 
angular momentum L = 0, and v) the mutual break-up channel with L = 2. Both of the 
channels iv) and v) are coupled to the channels ii) and iii) by the El operator of the target 
and of the projectile, respectively. The adiabatic potential AVad{R) is given as the lowest 
eigenvalue of the 5x5 coupling matrix at each R. Here we neglect the dipole-dipole term in 
the interaction, which we assume to be much smaller than the monopole-dipole term. 

Let us now numerically estimate the effect of the virtual break-up coupling on astro- 
physical fusion reactions. We first consider the effect of deuteron break-up on the d+^Re 
reaction. The break-up Q-value is e—2.22 MeV, and Eq.(6) leads to the total B{E1) f 
strength of 0.558 (e^ fm^). The dashed line in Fig. 1 shows the enhancement factor / of the 
penetrability, P{E)/ Po{E), as a function of the center of mass energy Ec.m., where Po{E) 
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is the penetrability of the bare Coulomb interaction Vo{R). We take Rq — 4.3 fm for the 
inner turning point [18]. We see that the enhancement factor slowly decreases as the energy 
decreases. The value of the enhancement factor is about 0.21 % at £^c.m. — 5.8 keV, which 
is smaller than the effect of vacuum polarization [18] by one order. Therefore, the effect 
of the dipole polarizability of deuteron seems to be negligible as compared to the vacuum 
polarization effect. The dotted line shows the effect of ^He break-up, where the binding 
energy is 5.49 MeV from the threshold of the d + p system. Although this effect is much 
smaller than the deuteron break-up effect, the combined effect of the mutual excitations 
increases the penetrability in a non-negligible way (see the solid line). Figure 2 shows the 
effect of dipole break-up of ^Li nucleus (into a-\-t) on the p-|-^Li reaction. For this system, 
the El effective charge is small, and the effect of break-up is even much smaller than the 
d+^Hc system. Notice that the El effective charge vanishes for a similar projectile, ^Li, 
which predominantly breaks into a + d. 

In summary, we have studied the effect of virtual Coulomb break-up process of colliding 
nuclei (i.e., the dipole polarizability) on astrophysical fusion reactions. For the deuteron 
break-up, we have found that the enhancement of the tunneling probability is about 0.2 
% for the d + '^He system. The effect is much smaller for the ^Li break-up in the p+^Li 
system, where the enhancement factor was found to be about 2.7x10"^ %. Therefore, the 
break-up effect alone does not resolve the large screening puzzle. We have a feeling that we 
have almost exhausted the list of small effects in astrophysical reactions. Of course, there 
are still some exotic effects such as the deformation of proton [38] or the color van der Waals 
force [39], but these effects should be extremely small in the astrophysical reaction. We 
may now be at a stage where the atomic physics based model has to be re-examined with a 
more careful and consistent treatment of few-body dynamics of charged particles including 
electrons. 

We thank the IPN Orsay for their warm hospitality where this work was carried out. 
K.H. also acknowledges support from the Kyoto University Foundation. This work was 
supported in part by the U.S. National Science Foundation Grants No. PHY-0244384 and 
INT-0070889 at the University of Wisconsin. 
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FIG. 1. The effect of virtual Coulomb break-up of colliding nuclei on the d+^He reaction. / 
is the enhancement factor of penetrability due to the break-up, measured from unity. The dashed 
and the dotted lines show the effect of break-up of the d and the "^Hc nuclei, respectively. The solid 
line is the combined effect of mutual break-up of both the projectile and the target nuclei. 



6 




7 



REFERENCES 



[1] H.J. Assenbaum, K. Langanke, and C. Rolfs, Z. Phys. A327, 461 (1987). 
[2] S. Engstler, A. Krauss, K. Neldner, C. Rolfs, U. Schroder, and K. Langanke, Phys. Lett. 
B202, 179 (1988). 

[3] F. Strieder, C. Rolfs, C. Spitaleri, and P. Corvisiero, Naturwissenschaften 88, 461 
(2001). 

[4] K. Langanke, in Advances in Nuclear Physics, edited by J.W. Negele and E. Vogt 

(Plenum, New York, 1993). 
[5] E.G. Adelberger et al, Rev. Mod. Phys. 70, 1265 (1998). 

[6] M. Aliotta, F. Raiola, G. Gyiirky, A. Formicola, R. Bonctti, C. Broggini, L. Campajola, 
P. Corvisiero, H. Costantini, A. D'Onofrio, Z. Fiilop, G. Gervino, L. Gialanella, A. 
Guglielmetti, C. Gustavino, G. Imbriani, M. Junker, P.G. Moroni, A. Ordine, P. Prati, 
V. Roca, D. Rogalla, C. Rolfs, M. Romano, F. Schiimann, E. Somorjai, O. Straniero, F. 
Strieder, F. Terrasi, H.P. Trautvetter, and S. Zavatarelli, Nucl. Phys. A690, 790 (2001). 

[7] H. Costantini, A. Formicola, M. Junker, R. Bonetti, C. Broggini, L. Campajola, P. 
Corvisiero, A. D'Onofrio, A. Fubini, G. Gervino, L. Gialanella, U. Greife, A. Gugliel- 
metti, C. Gustavino, G. Imbriani, A. Ordine, P.G. Prada Moroni, P. Prati, V. Roca, D. 
Rogalla, C. Rolfs, M. Romano, F. Schiimann, O. Straniero, F. Strieder, F. Terrasi, H.P. 
Trautvetter, and S. Zavatarelli, Phys. Lett. B482, 43 (2000). 

[8] G. Fiorentini, C. Rolfs, F.L. Villante, and B. Ricci, Phys. Rev. C67, 014603 (2003). 

[9] L. Bracci, G. Fiorentini, V.S. Melezhik, G. Mezzorani, and P. Quarati, Nucl. Phys. 
A513, 316 (1990). 

[10] N. Takigawa, K. Hagino, M. Abe, and A.B. Balantekin, Phys. Rev. C49, 2630 (1994). 

[11] N. Takigawa, K. Hagino, and M. Abe, Phys. Rev. C51, 187 (1995). 

[12] A.B. Balantekin and N. Takigawa, Rev. Mod. Phys. 70, 77 (1998). 

[13] T.D. Shoppa, S.E. Koonin, K. Langanke, and R. Seki, Phys. Rev. C48, 837 (1993). 

[14] F.C. Barker, Nucl. Phys. A707, 277 (2002). 

[15] C. Spitaleri, S. Typel, R.G. Pizzone, M. Ahotta, S. Blagus, M. Bogavac, S. Cherubini, 

P. Figuera, M. Lattuada, M. Milin, D. Miljanic, A. Musumarra, M.G. Pehegriti, D. 

Rendic, C. Rolfs, S. Romano, N. Soic, A. Tumino, H.H. Wolter, and M. Zadro, Phys. 

Rev. C63, 055801 (2001). 
[16] A. Musumarra, R.G. Pizzone, S. Blagus, M. Bogovac, P. Figuera, M. Lattuada, M. 

Mihn, D. Miljanic, M.G. Pellegriti, D. Rendic, C. Rolfs, N. Soic, C. Spitaleri, S. lypel, 

H.H. Wolter, and M. Zadro, Phys. Rev. C64, 068801 (2001). 
[17] A. Tumino, C. Spitaleri, A. Di Pietro, P. Figuera, M. Lattuada, A. Musumara, M.G. 

Pellegriti, R.G. Pizzone, S. Romano, C. Rolfs, S. Tudisco, and S. Typel, Phys. Rev. 

C67, 065803 (2003). 

[18] A.B. Balantekin, C.A. Bertulani, M.S. Hussein, Nucl. Phys. A627, 324 (1997). 

[19] K. Hagino and A.B. Balantekin, Phys. Rev. C66, 055801 (2002). 
[20] M.S. Hussein, Phys. Lett. 71B, 249 (1977). 
[21] M.S. Hussein, Phys. Rev. C19, 807 (1979). 

[22] Z.E. Switkowski, R.M. Wieland, and A. Winther, Phys. Rev. Lett. 33, 840 (1974). 
[23] D. Galetti and M.A. Candido Ribeiro, Phys. Rev. C50, 2136 (1994). 
[24] L.C. Chamon, D. Pereira, M.S. Hussein, M.A. Candido Ribeiro, and D. Galetti, Phys. 
Rev. Lett. 79, 5218 (1997). 

8 



[25] A.B. Balantekin, J.F. Beacom, and M.A. Candido Ribciro, J. Phys. G24, 2087 (1998). 

[26] G. Bliige and K. Langanke, Phys. Rev. C41, 1191 (1990). 

[27] C.H. Dasso, M. Lozano, and A. Vitturi, Phys. Rev. A44, 4743 (1991). 

[28] J.C. Fernandas, R. Crespo, F. Nunes, and I.J. Thompson, Phys. Rev. C59, 2865 (1999). 

[29] K. Ogata, M. Yahiro, Y. Iseri, and M. Kamimura, Phys. Rev. C67, 011602(R) (2003). 

[30] C.H. Dasso and A. Vitturi, Phys. Rev. C50, R12 (1994). 

[31] K. Hagino, A. Vitturi, C.H. Dasso, and S.M. Lenzi, Phys. Rev. C61, 037602 (2000). 
[32] A. Diaz-Torres and I.J. Thompson, Phys. Rev. C65, 024606 (2002). 
[33] M.S. Hussein, M.P. Pato, L.F. Canto, and R. Donangelo, Phys. Rev. C47, 2398 (1993). 
[34] W.H.Z. Cardenas, L.F. Canto, R. Donangelo, M.S. Hussein, J. Lubian, and A. Ro- 

manelh, Nud. Phys. A703, 633 (2002). 
[35] C.A. Bertulani, G. Baur, and M.S. Hussein, Nud. Phys. A526, 751 (1991). 
[36] K. Alder, A. Bohr, T. Huus, B.R. Mottelson, and A. Winther, Rev. Mod. Phys. 28, 432 

(1956). 

[37] M.S. Hussein, V.L.M. Franzin, R. Franzin, and A.J. Baltz, Phys. Rev. C30, 184 (1984). 
[38] Sec e.g., G.A. Miller, e-print: nud-th/0304076. 

[39] M.S. Hussdn, C.L. Lima, M. Pato, and C.A. Bertulani, Phys. Rev. Lett. 65, 839 (1990). 



9 



